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Atastiact 

The  cause  of  a  discrepancy  in  the  CCIR  Report  322-3  radio  noise  model  is  reported.  The  basis  for  this  discrepancy  resulis  from 
the  procedure  used  to  prepare  the  measured  noise  data  for  the  determination  of  a  global  numerical  representation  of  the  I  Mh/  data 
The  development  process  of  CCIR  Report  322-3  is  presented  so  that  the  cause  of  the  error  can  be  more  fully  understood  Then  the 
cause  itself  of  the  error  in  the  I  Mhz  noise  model  is  discussed.  Next  gei^raphical,  frequency  dependence,  and  interpolation  effects 
are  reviewed.  Finally,  brief  recommendations  of  a  course  of  action  for  the  development  of  an  iniproved  1  Mhz  model  are  iiivcn 


Introduction 

Recently,  Bowen  and  Fraser-Smith  (1992)  made  a  comparison  of  measured  32  Khz  radio  noise  amplitudes  with  the  CCIR  Report 
322-3  (1988)  noise  model  predictions.  They  found  that  the  greatest  discrepancies  between  the  measured  and  predicted  amplitudes 
were  observed  at  the  two  northern  high  latitude  stations  (Sandre  StramQord  and  Thule.  Greenland),  where  on  some  occasions  the 
predicted  values  were  nearly  five  times  greater  than  the  measured  values.  There  was  moderately  good  agreement  between  the 
measured  and  predicted  values  at  a  southern  high  latitude  site  (Arrival  Heights,  Antarctica).  The  best  agreement  was  observed  at 
middle  to  low  latitudes.  The  data  used  to  make  these  comparisons  were  measured  by  a  ELF/VLF  measutcraent  system  (Fraser-Smith 
et  al..  1987). 

CCIR  Report  322-3  is  an  output  document  of  the  CCIR  XVlth  Plenary  Assembly  held  in  Dubrovnik,  Yugoslavia  in  1986.  It 
was  produced  by  the  International  Working  Party  (fWP)  6/2  of  CCIR  Study  Group  6  using  the  work  of  Spaulding  and  Washburn 
(1985). 

Because  of  the  wide  acceptance  of  the  CCIR  Report  322-3  radio  noise  model,  an  attempt  was  made  to  determine  the  cause  of 
this  discrepancy.  This  error  resulted  from  the  procedure  used  to  prepare  the  measured  noise  data  for  the  determination  of  a  global 
numerical  representation  of  the  I  Mhz  data.  Since  the  contour  plots  of  the  I  Mhz  radio  noise  in  CCIR  Report  322-3  were  in  turn 
generated  from  the  numerical  representation  thus  developed,  they  are  also  in  error. 

The  Process  Used  To  Develop  TIic  CCIR  Report  322-3  Radio  Noise  Model 

CCIR  Report  322  (1964)  was  developed  using  data  available  through  October  1961.  Data  from  the  original  worldwide  network 
of  recording  stations  continued  to  be  measured  through  1966.  Many  years  of  data  from  10  Soviet  measurement  locations  became 
available  along  with  data  from  Tliailand  for  March  1966  through  February  1968  (Chindahpom  and  Younker,  1968).  All  of  this  data 
was  analyzed  and  updated  set  of  atmospheric  radio  noise  estimates  produced,  essentially  in  the  CCIR  Report  322  format  (Spaulding 
and  Washburn,  1985). 

The  discussion  here  of  the  development  of  the  new  CCIR  Repon  322-3  radio  noise  model  is  based  on  the  work  of  Spaulding 
and  Washburn  (1985). 

The  New  Data 

The  original  worldwide  network  locations  and  new  locations  are  given  in  Table  1  (Table  1  in  Spaulding  and  Washburn  (1985)). 

In  the  development  of  the  new  noise  model,  data  from  Thule,  Greenland  and  Byrd  Station,  Antarctica  were  not  used.  This  data 
was  not  used  because  it  was  assumed  that  it  was  generally  contaminated  by  high  levels  of  man-made  noise. 

For  a  number  of  years  the  Soviet  Union  operated  a  network  of  ten  noise  measurement  stations.  Raw  data  were  available  on 
microfilm  from  the  World  Data  Center  (National  Oceanic  and  Atmospheric  Administration,  Boulder.  CO  80303)  for  periods  of  time 
from  mid-1958  through  1965,  coincident  with  the  measurements  of  the  worldwide  network.  The  parameters  that  were  measured 
were  different  from  those  discussed  for  the  worldwide  network  above.  The  conversion  of  these  parameters  to  parameters  equivalent 
to  that  of  the  worldwide  network  are  discussed  by  Spaulding  and  Washburn  (1985).  The  measurement  frequencies  and  other 
information  for  each  of  the  Soviet  measurement  locations  are  also  summarized  in  Spaulding  and  Washburn  (1985). 


Table  1.  Atmospheric  noise  measurement  locations 


Balboa,  Canal  Zone 

79. 5  W 

9.0N 

Alma  Ata.USSR 

76.92E 

43.25N 

Bill,  Wyoming 

105i2W 

43. 2N 

Ashkhabad,  USSR 

58  .U: 

37  92N 

Boulder.  Colorado 

105. IW 

40.  IN 

Irkutsk,  USSR 

104  5  F, 

52  ON 

Byrd  Station,  Antarctica 

I20.W 

80.0S 

Khabarovsk.  USSR 

135  01: 

50  ON 

Cook,  Australia 

130.4E 

30,6S 

Kiev.  USSR 

.30  3  E 

.50  72N 

Enkoping,  Sweden 

17.3E 

59.5N 

Laem  Chabang, 

100  9E 

13  05N 

Front  Royal,  Virginia 

7g,2W 

38.8N 

Moscow,  USSR 

37,32E 

55.47N 

Ibadan,  Nigeria 

3.9E 

7.4N 

Murmansk,  USSR 

35.0E 

69.0N 

Kekaha,  Hawaii 

I59.7W 

22.0N 

Simferopol,  USSR 

34,03E 

45.02N 

New  Delhi,  India 

77.3E 

28.8N 

Sverdlovsk.  USSR 

6I.07E 

56.73N 

Ohira,  Japan 

140.5E 

35. 6N 

Tbilisi.  USSR 

40.0E 

4I,72N 

Pretoria,  South  Africa 

28.3E 

25. 8S 

Rabat,  Morocco 

6.8W 

33.9N 

San  Jose,  Brazil 

45. 8W 

23.3S 

Singapore 

103. 8E 

1.3N 

Thule,  Greenland 

68. 7W 

76. 6N 

Analysis  of  Ihe  Soviet  Data 

The  analysis  involved  determining,  at  each  frequency,  for  each  3-nionth  period  and  4-hour  time  block,  the  median  value  of  all 
the  data.  These  median  values  at  the  various  frequencies  were  then  used  to  determine  the  approximate  I  Mhz  value.  This  value 
was  then  used  to  obtain  a  correction  value  to  the  CCIR  Report  322  value.  Figure  I  (Figure  8  in  Spaulding  and  Washburn  (1985)) 
shows  an  example  for  Moscow  for  June,  July,  August  (Northern  Hemisphere  Summer)  and  1600-2000  hours.  Note  that  the  noise 
curve  for  a  quiet  receiving  location  falls  considerably  below  the  analyzed  median  values.  A  computer  algorithm  was  developed  that 
determined  the  atmospheric  noise  frequency  variation  curve  that  "best"  fit  the  data.  However,  since  the  median  value  at  some 
frequencies  was  based  on  much  more  data  than  the  value  of  the  oilier  frequencies  for  a  location,  time  block,  and  season  (due  to 
missing  data  and  some  frequencies  being  stressed  at  some  locations);  this  'fitting*  process  was  generally  done  by  hand  (visually). 
On  Figure  1,  the  'best'  fitting  frequency  law  curve  was  determined  to  be  72  Db.  The  CCIR  Report  322  value  is  65  dB.  Hence, 
a  value  of  +1  dB  can  be  used  to  correct  the  atmospheric  noise  predicted  by  CCIR  Report  322  at  Moscow  during  the  summer  at 
1600-2000  hours.  Figure  2  (Figure  9  in  Spaulding  and  Washburn  (1985))  shows  an  example  for  Moscow  for  the  period  November, 
December,  January,  0800-1200  hours.  Atmospheric  noise  would  be  expected  to  be  low  during  this  period  (winter  morning)  and 
could  possibly  be  contaminated  by  man-made  noise  at  the  higher  frequencies.  For  the  higher  frequencies  350  kHz  and  above,  the 
Figure  shows  a  typical  man-made  noise 
curve  at  a  level  to  be  expected  for  a  quiet 
receiving  site.  Because  of  this  contamina¬ 
tion  possibility,  the  lower  frequencies  were 
used  to  determine  the  frequency  law  curve. 

In  this  cast  the  frequency  law  curve  for  31 
dB  was  determined.  The  CCIR  Repon  322 
value  was  29  dB,  resulting  in  a  required 
correction  of  -t-2  dB. 

Corrections  to  CCIR  Report  322  I  MHz 
Values 

The  procedure  illustrated  above  for 
Figures  I  and  2  for  determining  the  cor¬ 
rections  to  be  made  to  CCIR  Report  322 
was  followed  by  Spaulding  and  Washburn 
(1985)  to  obtain  corrections  for  each  loca¬ 
tion  and  for  each  3-monlh/4-hour  lime 
block.  Tables  2  gives  the  corrections  deter¬ 
mined  using  this  procedure  for  the  Decem¬ 
ber,  January,  February  season.  Similar 
tables  for  the  other  three  seasons  are  given 
in  Spaulding  and  Washburn  (1985).  The 
table  contains  a  correction  for  time  block 
and  for  each  station  in  Table  I  except  for 
those  noted  below.  The  "correction"  is  the 
difference  between  Ihe  CCIR  Report  322  I 
MHz  value  and  the  corresponding 
value  determined  using  the  above  procc-  f'e“re  1.  Determination  of  1  MHz  F„  value  for  Moscow,  June,  July,  August,  1600- 
dure  from  the  data.  (Spaulding  and  Washburn,  1985) 

For  certain  stations  listed  in  Table  1 ,  no  corrections  were  determined.  No  correction  values  were  obtained  for  Thule,  Greenland 
and  Byrd  Station,  Antarctica  because  it  was  assumed  that  the  data  v-s  contaminated  by  man-made  noise.  Since  there  verr  no  data 
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trom  Ibadan,  Nigeria  past  the  data  used  in  the  development  of  CCIR  Report  322;  no  correction  was  used  loi  Ihadoii  l-ksause  ilie 
corrections  for  Bill,  Wyoming  and  Boulder,  Colorado  were  essentially  identical,  only  corrections  loi  Boiildei  weie  used 
Corrections  were  used  for  only  6  Soviet 
locations  rather  than  10.  Simferopol, 

Sverdlovsk,  Tbilisi,  and  Kiev  ehher  had 
only  small  amounts  of  usable  low  frequen¬ 
cy  data  necessary  lo  determine  the  proper 
I  MHz  value  oi  were  close  to  other 
measurement  locations.  The  data  at  these 
four  Soviet  locations  were  analyzed  lo 
ascertain  that  the  corrections  agreed  with 
those  used  at  nearby  locations,  namely 
Moscow  and  Ashkhabad.  For  Murmansk; 
for  the  March,  April,  May  season;  and  for 
the  four  time  blocks  0400-0800,  0800- 
1200,  1200-1600,  and  1600-2000  hours, 
the  dauj  were  highly  irregular  and  confus¬ 
ing.  Hence,  Spaulding  and  Washburn 
(1985)  decided  not  to  attempt  lo  obtain  any 
correction  values  for  Murmansk  for  these 
four  periods. 


Because  the  original  contour  maps  in 
CCIR  Report  322  were  produced  directly 
from  a  grid  of  equally  spaced  84  longitude 
by  100  latitude  points  for  each  time 
block/season,  the  next  step  in  the  analysis 
was  to  translate  the  correction  data  for  the 
nineteen  sites  to  this  same  84  by  100  lat¬ 
tice.  Then  data  for  new  noise  model  was 
found  simply  by  adding  the  corrections 
point-by-point  for  the  grid  point  of  original 
data. 

To  do  this  Spaulding  and  Washburn  (1985)  utilized  an  interpolation  method  due  to  Dr.  Charles  L.  Lawson  (1982,  1984)  for 
interpolating  scattered  data  over  a  sphere.  This  method  first  constructs  a  triangular  grid  over  a  surface  (in  this  case  the  Earth)  using 
a  given  set  of  points  as  vertices  (in  this  case  the  19  data  points  in  Tables  2  and  Tables  3  through  S  in  Spaulding  and  Washburn 
(I98S) .}.  Then  continuous  first  partial  derivatives  are  estimated  by  the  method  at  each  vertex  using  local  quadratic  least  squares 
approximations  to  given  data  values  at  nearby  vertices.  The  method  for  interpolation  then  uses  six  Hermite  cubic  interpolations 
along  arcs  of  great  circles.  Figure  3  is  an  example  of  the  resultant  contour  maps  of  the  24  (four  3-month  periods,  six  4-hour  time 
blocks)  100  by  84  correction  grids  produced  by  Spaulding  and  Washburn  for  December,  January,  February  for  12001600  hours 
(The  other  23  maps  can  be  found  in  Spaulding  and  Washburn  (1985)  as  Tables  10  to  12  and  as  Tables  14  to  33).  These  maps  show 
the  changes  to  the  CCIR  Report  322  lo  be  made.  There  are  substantial  corrections  in  some  areas  as  could  also  be  seen  by  examining 
•Tables  2  and  Tables  3  through  5  in  Spaulding  and  Washburn  (1985).  The  correction  maps  are  presented  m  terms  of  three-month 
periods  rather  than  as  season  (which  results  in  a  discontinuity  at  the  equator)  as  in  CCIR  Report  322. 

The  New  1  MHz  Fj,^  Values 

The  new  I  MHz  Fj^  data  values  for  constructing  a  new  1  MHz  noise  model  where  obtained  by  Spaulding  and  Washburn  (1985) 
by  adding  each  of  the  84  longitude  by  100  latitude  grid  of  correction  values  lo  the  corresponding  original  84  longitude  by  100 
latitude  grid  data  values  from  which  CCIR  Report  322  was  constructed. 

In  developing  a  numerical  representation  for  these  new  1  MHz  Fj,„  maps,  Spaulding  and  Washburn  used  the  method  used  by 
Lucas  and  Harper  (1965).  Spaulding  and  Washburn  give  details  on  how  these  coefficients  were  obtained. 

These  numerical  maps  represent  a  "smoothed"  version  of  the  original  data  and  are  the  new  I  MHz  Fjj„  worldwide  atmospheric 
noise  estimates.  Figures  34-57  in  Spaulding  and  Washburn  (1985))  are  contour  plots  of  these  estimates.  Upon  comparing  the 
numerical  representation  thus  obtained  for  each  of  the  8400  original  data  points  (84  x  l(X)  grids)  for  each  of  the  24  numerical  maps. 
Spaulding  and  Washburn  found  an  rms  variation  that  ranged  from  0.88  dB  to  2.37  dB  over  the  24  maps  with  an  average  rms 
variation  of  1 .52  dB  and  with  a  maximum  deviation  of  6.7  dB  (all  maps  considered,  i.e.,  24  x  8400  points).  The  contour  plots  in 
CCIR  Report  322-3  are  similar,  except  the  CCIR  IWP  6/2  gathered  the  plots  together  by  season  rather  than  by  months.  This 
resulted  in  the  discontinuity  at  the  equator  shown  in  CCIR  Report  322-3. 

The  Probable  Source  of  the  Error  in  the  CCIR  Radio  Noise  Model 

The  probable  source  of  the  error  in  the  CCIR  radio  noise  model  is  most  likely  due  to  the  non-use  of  correction  factors  for  Thule, 
Greenland;  Byrd  Station,  Antarctica;  Ibadan,  Nigeria;  and  Bill,  Wyoming.  In  the  case  of  Thule  and  Byrd  Station,  data  from  these 
sites  were  used  in  the  original  CCIR  noise  model.  One  would  anticipate  that  if  the  data  were  contaminated  bv  man-made  noise  that 
a  negative  correction  factor  would  be  th"*  Mkelv  case  Thus,  it  would  seem  dcsirabl'-  to  deli-iuine  cc  rertion  factors  for  tiicse  sites 
as  was  done  for  the  other  sites.  Tliat  is,  the  data  at  lower  frequencies  could  have  been  used  to  obtain  correction  factors  at  I  MHz. 


Figure  2.  Determination  of  1  MHz  F„  value  for  Moscow,  December,  January, 
February,  1000-1200  hours  (Spaulding  and  Washburn,  1985) 
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CertaiTily,  there  Is  evidence  that  the  measureineius  at  Thule  at  2.5  and  5.0  MHr  were  contaminated  by  manmade  noise  (Ikiinan. 
1962).  However,  Herman  (1963)  showed  that  Byrd  Station  is  an  exceptionally  kjuiet  location  and  would  be  a  good  sue  loi  iiiaking 
a  variety  of  radio  measurements  requiring  a  low  noise  background.  At  Thule  man-made  noise  on  2  5  and  5  0  MHz  appeared  to  be 
57  and  49  dB  above  kTB,  respectively,  while  at  Byrd  Station  the  values  were  about  20  and  12  dB,  respectively  (Herman,  I9t>4) 
These  values  were  estimated  from  data  taken  during  a  PCA,  when  atmospheric  noise  was  absent.  The  corresponding  values  for 
a  quiet  rural  site  at  2.5  and  5.0  MHz  are  43  and  34  dB,  respectively.  Both  sites  are  affected  by  galactic  noise  at  10  0  and  20  0 
MHz.  In  the  case  of  Ibadan,  no  addition  data  was  available  beyond  the  last  date  for  data  used  in  the  original  CCIR  322  noise  model 
Not  even  a  correction  factor  of  zero  was  used  to  mainuin  the  status  quo  was  used  for  these  three  sites.  In  the  case  of  Bill,  data 
was  not  included  because  the  correction  factors  obtained  for  it  nearly  the  same  as  for  Boulder,  Colorado,  which  was  in  close 
proximity  to  Bill,  Wyoming. 

Table  2.  Corrections  (dB)  to  CCIR  Report  322  I  MHz  values  for  December,  January,  and  February  (Spaulding  and 
Washburn  (1985)) 


Figure  3.  Corrections  (dB)  to  original  CCIR  Report  322  1  MHz  estimates,  December,  January,  February,  12(X)-16(X)  hours 
(Spaulding  and  Washburn,  1985) 


The  nel  result  of  not  including  coriection  factors  (not  even  zero  to  maintain  the  stains  quo)  lor  these  limi  locations  was  that 
the  inteipolalion  algonihm  used  to  deternilne  the  100  latitude  by  84  longitude  data  points  supplied  crioiuous  correction  las  lois  at 
these  sites.  Tables  3  through  6  give  the  eiror  In  the  correction  factor  for  these  four  sites.  Tor  Thule.  Byrd  Station,  and  Ibadan,  the 
erroi  is  the  difference  between  the  correction  factors  given  in  Tigures  10  through  33  in  Spaulding  and  Washburn  and  zcio  lor  tin- 
status  quo.  Tor  Bill  the  error  is  the  difference  between  those  given  in  these  figures  and  the  correction  lacioi  input  tot  Houklei. 
Colorado.  It  was  anticipated  th.at  the  errors  at  3'hule  and  Byrd  station  might  be  large,  but  it  was  a  suipnse  see  the  magnitude  oi 
the  errors  at  Ibadan.  Tor  Thule  the  maximum  and  minimum  errors  in  the  correction  contours  from  zero  siaius  quo  weie  1(1  1  and 
-10.8  dB,  respectively.  For  Ibadan  the  niaximum  and  minimum  errors  were  1 2.5  and  -1.5  dB.  res|>eciivelv  I  or  Byrd  Station  die 
maximum  and  minimum  errors  were  12.0  and  .3  dB,  respectively.  The  error  for  Bill,  Wyoming  is  wuhin  ilic  nns  crior  oi  the 
luimcrical  maps  of  CCIR  Report  322  3 

Table  3.  Interpolation  errors  (dB)  for  select  measurement  locations  for  December,  January,  February 


LOCATION 

LOCAL  TIME  1 

NAME 

112011 

08-12 

t  12-16 

16-20 

20-24  1 

Thule 

7  0 

3  5 

6.0 

Ibadan 

-0  4 

IggHI 

5.5 

3  0 

-1  5 

Bill 

-0  8 

-0.3 

-I.O 

-1.0  1 

Byrd  Station 

3.0 

0.0 

5  5 

6  0 

l5  1 

Table  4.  Interpolation  errors  (dB)  for  select  measurement  locations  for  March,  April,  and  May 


LOCATION 

LOCAL  TIME  j 

NAME 

00-04 

[  04-08 

08-12 

12-16 

16-20 

20-24  1 

Thule 

3,0 

-3.0 

-4.1 

3.0 

1  Ibadan 

3.0 

lyw 

8.2 

0.8 

5.0 

6.0 

Bill 

-1.0 

-1.0 

-6.6 

j  Byrd  Station 

4.0 

8.0 

12.0 

7.3 

10.2 

7.6  1 

Table  5.  Interpolation  errors  (dB)  for  select  measurement  Irxations  for  June,  July,  and  August 


LOCATION 

LOCAL  TIME  J 

NAME 

00-04 

04-08 

08-12  1  12-16 

1  16-20 

1  20-24  [ 

Thule 

■1.3 

6.0 

1.5  i  6.0 

HjXHi 

Ibadan 

5.3 

7.4 

10.6  1  0.0 

4.5 

■qH' 

Bill 

-0,5 

-1.0 

-1.0  i  -1.0 

-1.0 

i  -0.7  1 

Byrd  Station 

8.3 

8.0 

12.0  i  7.3 

10.2 

^  ’0 

Table  6.  Interpolation  errors  (dB)  for  select  measurement  locations  for  September,  October,  November 


1  LOCATION 

LOCAL  TIME  | 

1  NAME 

00-04 

j  04-08 

1  08-12 

1  12-16 

16-20 

20-24 

1  Thule 

-6.0 

3.0 

I  0.0 

1  3.0 

3.0 

-2.5 

1  Ibadan 

5.5 

10.7 

!  12-5 

1  5.2 

8.8 

5.3 

Bill 

-1.4 

-1,0 

j  -1.0 

j  -1.0 

-0.7 

1  Byrd  Station 

3.0 

4.9 

1  12.0 

i  9.0 

3.0 

3.0  1 

Geographical  Effects 

To  see  the  geographical  effect  of  these  errors,  contour  plots  were  made  of  the  errors  for  each  time  block  and  season  For  tlic 
19  locations  used  to  determine  the  correction  factor  contours  in  Figures  10  through  33  in  Spaulding  and  Washburn  (1985),  the  error 
was  assumed  to  be  zero.  It  was  assumed  that  the  interpolation  process  gave  the  same  values  for  the  correction  factors  at  these  sites 
as  was  input.  The  values  given  in  Tables  3  through  6  were  used  for  Ute  other  four  sites.  The  graphics  program  Axum  was  used 
to  plot  the  errors  for  the  irregularly  spaced  data  points.  Axum  rarely  was  able  to  determine  a  grid  spacing  greater  than  fifty  by  fifty 
for  its  internal  interpolation  for  the  contour  plotting.  Examples  of  Ihe  results  of  Sailors  (1993)  are  presented  in  Figures  4  through 
7.  Note  that  the  longitudes  are  positive  for  degrees  East  of  Greenwich  (zero  degree)  and  negative  for  degrees  West  of  zero  degree. 
Positive  latitudes  are  Northern  latitudes,  and  negative  latitudes  arc  Southern  latitudes.  The  locations  of  the  four  sites  for  which  no 
correction  factors  were  used  in  the  interpolation  are  approximately  given  in  the  figures.  Examination  of  these  figures  reveals  Ihai 
the  geographical  extent  of  the  error  is  not  confined  to  the  measurement  location  but  in  fact  in  some  cases  is  very  large.  This  is 
especially  true  in  the  Northern  and  Southern  high  latitudes,  the  Arabian  Peninsula,  Northern  Africa,  and  the  Mid-Atlantic  areas. 
Note  that  the  geographical  extent  of  the  error  is  both  seasonally  and  diumaJly  dependent. 
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Frequency  DeiJendeiice  Effects 

As  the  CCIR  Report  322-3  atmospheric  noise  model  can  be  used  at  frequencies  from  10  kHz  through  30  MHz,  ii  is  ini|xin.iiii 
to  be  able  to  translate  an  error  m  the  I  MHz  model  to  an  error  at  any  frequency  in  this  frequency  range  lo  determine  tins 
relationship,  the  frequency  dependence  model  was  used,  assuming  no  error  inherit  in  it  This  model  obtains  the  noise  ai  .iny 
frequency  by  putting  in  the  1  MHz  value  for  a  particular  time  block  and  season  into  a  set  of  curves  parametric  in  the  I  MHz  saUic 
For  each  time  block  and  seasons,  there  is  a  range  of  possible  input  I  MHz  values  for  which  the  atmospheric  noise  can  be  obtained 
depending  on  the  location  of  the  receive  site.  To  obtain  the  error  at  an  arbitrary  frequency,  time  block,  and  season;  each  ol  the 
parametric  values  were  penurbed  by  10  and  -10  dB  errors  and  the  resultant  errors,  respectively,  were  deierinmed  Hie  si.ii,siic.il 
average,  high  and  low  values  were  then  obtained.  It  was  found  that  the  error  at  each  of  the  35  frequencies  Irom  10  kHz  iliumgli 
30  MHz  used  did  no.  depend  significantly  on  the  parametne  curve  value.  However,  the  error  was  diurnally  dependent  T  he 
maximum  errors  for  each  frequency  occurred  during  the  local  lime  daytime,  and  the  minimum  error  errors  occurred  during  the 
nighttime.  Figure  8  shows  the  average, 

high,  and  low  error  as  a  function  of  r - ■ 

frequency  and  I  MHz  error.  The  figure 
shows  the  error  for  - 10  dB  to  be  a  mirror 
image  of  that  for  a  10  dB  error. 


Interpolation  Effects 


The  accuracy  of  the  interpolation  itself 
is  affected  by  the  lack  of  inclusion  of 
data  from  the  four  measurement  loca¬ 
tions.  First,  the  number  of  triangles  and 
number  of  edges  of  these  triangles  is 
reduced.  Lawson  (1977)  gives  equations 
for  both  the  number  of  triangles  n,  and 
number  of  edges  n,.  as  a  function  the 
number  distinct  points  n  in  the  set  S,  the 
number  of  points  n^,  in  S  on  the  boundary 
of  the  convex  hull  of  S,  and  the  numb^ 
of  points  nj  in  the  interior  of  the  convex 

hull  of  S  so  that  n  =  n|j  -F  nj.  In  the _ 

case  of  data  points  on  a  sphere,  nj,  =  0  Figure  4.  Geographical  I  MHz  atmospheric  noise  model  error,  December,  January, 
(i.e.,  there  are  no  boundary  points).  For  February,  1200-1600  hours 
n  =  19,  n,  =  36  and  n^  =  54.  For  n  •= 

23,  n,  =  44  and  n,.  =  66.  Thus  there  is 
an  122.22  percent  increase  in  triangle 
edges  by  adding  the  four  points.  Second, 
the  number  of  input  data  points  affects 
how  the  c'  surface  interpolation  deter¬ 
mines  the  triangles  themselves.  Each  pair 
of  triangles  forms  a  quadrilateral.  The 
Lawson  criterion  requires  that  the  small¬ 
est  of  the  six  angles  in  the  two  triangles 
of  a  quadrilateral  should  be  larger  for 
this  division  of  a  convex  quadrilateral 
than  that  given  by  the  other  diagonal 
(Ripley,  1981).  Thus  when  the  four  data 
measurement  locations  were  not  used  in 
the  model  development,  the  triangulation 
was  considerably  changed  and  may  not 
have  been  optimum.  Further,  Akima 
(1984)  has  shown  that  poor  estimates  of 
partial  derivatives  usually  occurs  when  a 
thin  (or  slim)  triangle  is  involved  in  the 
interpolation,  which  affects  the  accuracy 
of  the  interpolation.  This  case  was  more 
likely  to  have  occurred  when  the  four 
data  locations,  particularly  the  two  high  latitude  sites,  were  left  out. 

Several  of  the  contour  plots  giving  the  correction  factors  in  Spaulding  and  Washburn's  Figures  10  through  33  have  values  larger 
or  smaller  than  the  extreme  values  used  in  the  interpolation  to  obtain  these  figures.  The  figures  with  the  minimum  extreme  less  than 
the  input  minimum  extreme  usually  had  this  extreme  occurring  at  one  or  two  of  the  same  locations.  These  two  locations  were  cither 
in  the  very  most  Easiem  pan  of  Asia  or  in  the  Nonh  Atlantic  just  south  of  Greenland.  These  last  two  extremes  likely  would  not 
have  occurred  had  the  data  at  Thule.  Greenland  been  used;  in  this  case  the  triangulation  would  have  been  considerably  different. 
The  figures  with  the  maximum  extreme  greater  than  input  maximum  extreme  usually  had  this  extreme  occurring  at  one  or  two  of 
the  same  locations.  The  first  location  was  usually  in  the  vicinity  of  Guam.  An  example  is  shown  in  Figure  3  above.  The  second 
location  was  in  the  South  Atlantic  between  South  America  and  Africa.  This  latter  location  would  not  have  been  affected  by  the 
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Figure  5.  Geographical  I  MHz  atmospheric  noise  model  error,  March,  April,  May, 
16(X)-2C)00  hours 
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iriangulauon  if  (he  data  trom  Ibadan,  Nigeria  and  Byrd  Slalion,  Antarctica  had  been  included  Huu-ecei ,  ilie  eciieine  cKc.iiiiii;:  i:. 
the  vicinity  of  Guam  is  due  (he  gradients  occurring  in  the  input  data  in  Eastern  Asia  The  data  diere  g.ies  lumi  laigc  uej  jn^t 
corrections  in  the  North  to  Large  positive 
corrections  at  Singapore  The  interpola 
tion  program  extrapolated  to  obtain  the 
large  positive  correction  factors  in  the 
vicinity  of  Guam  Eurther,  in  this  case 
the  triangulatioii  would  not  have  changed 
if  (he  data  for  the  other  four  locations 
had  been  included  in  the  model  develop 
ment.  Whether  this  correction  is  valid  or 
not  is  unclear.  Addition  input  data  in  the 
vicinity  of  Guam  would  have  been  useful. 

Conclusions  and  Recommendations 

This  paper  has  presented  the  probable 
cause  lor  a  discrepancy  in  the  CCIR 
Report  322-3  radio  noise  model.  This 
report  has  not  determined  whether  there 
is  a  discrepancy  between  the  new  CCIR 
mrxlel  and  the  data  values  used  to  devel¬ 
op  it.  Nor  has  there  been  an  attempt  to 
determine  the  validity  of  the  measured 
data  values  used  to  develop  (he  model  or 
to  validate  the  model  against  any  other 
data. 

The  basis  for  (his  discrepancy  was  found  to  be  in  the  procedure  used  to  prepare  the  measured  noise  data  lor  the  deteiniimiiion 
of  a  global  numerical  representation  of  the  I  MHz  data.  The  procedure  followed  in  the  development  of  the  model  was  to  dcicrmme 
correction  factors  to  the  old  CCIR  model  for  each  measurement  site,  to  interpolate  these  corrections  to  a  100  latitude  by  84 
longitude  grid  for  each  time  block/season,  to  add  the  correction  factois  at  each  grid  point  to  corresponding  values  for  the  old  CCIR 
model,  and  finally  to  numerically  map  (he  resulting  data  for  each  time  block  and  season.  Nineteen  locations  were  used  in  the  final 
model.  Four  sites  used  in  the  original  CCIR  model  were  not  used.  These  include  Bill,  Wyoming;  Byrd  Station,  Antarctica;  Ibadan, 
Nigeria;  and  Thule,  Greenland.  As  no  correction  factors  were  obtained  for  these  locations  or  a  correction  factor  of  zero  used,  the 
interpolation  algorithm  used  to  obtain  the  100  latitude  by  84  longitude  grid  of  correction  factors  supplied  erroneous  vaiues.  For 
Bill,  Wyoming  the  result  is  not  too  serious;  but  for  the  other  three  sites,  the  error  is  at  some  seasons  and  time  of  day  serious.  For 
Thule,  Gre^land  the  maximum  and  minimum  errors  in  the  correction  contours  were  10. 1  and  -10.8  dB,  respectively.  For  Ibadan. 
Nigeria  the  maximum  and  minimum  errors  were  I2.S  and  -l.S  dB,  respectively.  For  Byrd  Station,  Antarctica  the  maximum  and 
minimum  errors  were  12.0  and  3.0  dB,  respectively.  Examination  of  the  geographical  extent  of  these  errors  reveals  that  the  error 
is  not  conftned  to  the  measurement  Icxation  but  in  fact  is  very  large.  (I  was  found  that  the  error  as  a  function  of  frequency  was 
diumally  dependent.  An  error  of  10  dB  at  1  MHz  was  more  serious  at  another  frequency  during  local  daytime  than  at  night  Finally, 
the  absence  of  tlie  data  locations  affected  the  accuracy  of  the  interpolation  itself. 

Because  of  the  errors  in  the  CCIR  Report  322-3  atmospheric  noise  model,  it  is  recommended  that  it  be  used  with  caution.  It 
is  most  accurate  in  Europe,  Asia,  the  Indian  Ocean,  the  Western  Pacific  from  Asia  to  the  date  line,  and  Australia.  It  is  most 
inaccurate  in  both  the  Northern  and  Southern  high  latitudes,  the  Arabian  Peninsula,  Northern  Afnca.  and  the  Mid-Atlantic  Ocean 
area.  For  applications  in  latter  list  of  areas,  (he  user  should  consider  using  the  original  CCIR  Report  372  model. 

A  three  step  process  should  be  followed  to  develop  a  new  1  MHz  atmospheric  noise  model.  First,  obtain  correction  factors  for 
additional  locations  to  increase  the  accuracy  of  the  interpolation.  Second,  test  the  method  of  interpolation  against  a  bench  mark. 
Third,  use  the  Zacharisen  and  Jones  (1970)  numerical  mapping  technique  applied  in  local  lime  to  develop  the  final  model.  Consider 
using  a  latitude  transformation  to  increase  the  accuracy  of  the  numerical  mapping  technique. 
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Figi..e  7.  Geographical  I  MHz  atmospheric  noise  model  error,June,  July, 
August,  0800-1200  hours 
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